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We present results from Monte Carlo simulations of the three-dimensional Gross-Neveu Model with a U(l) 
chiral symmetry at nonzero temperature. We provide evidence that the model undergoes a Berezinskii-Kosterlitz- 
Thouless transition in accordance with the dimensional reduction scenario. We also identify a regime in the high 
temperature phase in which the fermions acquire nonzero dynamical mass, analogous to the pseudogap behaviour 
observed in cuprate superconductors. 



1. Introduction 

The three-dimensional Gross-Neveu model 
(GNM3) has been proved to be an interesting and 
tractable model to study chiral phase transitions 
both numerically by means of lattice simulations 
and analytically in the form of large- Nf expan- 
sions. In this report we present results of numer- 
ical simulations of the U(Nf )v/-invariant GNM3 
with a C/(l) chiral symmetry at non-zero temper- 
ature. This model is described by the following 
continuum Euclidean Lagrangian density: 

£ = y i ($ + <j + i l5 n)y 2 + ^(a 2 +7T 2 ). (1) 

We treat ^>i, 'I'i as four-component Dirac spinors 
and the index i runs over Nf fermion species. 
The model is renormalizable in the 1/77/ expan- 
sion unlike in the loop expansion |Q. For suffi- 
ciently strong coupling at T — it exhibits sponta- 
neous chiral symmetry breaking and if we choose 
(ir) = the pion field tt is the associated Gold- 
stone boson. 

At leading order in l/Nf this model undergoes 
a second order chiral phase transition This 
conclusion is expected to be valid only when Nf 
is strictly infinite, i.e. when the fluctuations of the 
bosonic fields are neglected, otherwise it runs foul 
of the Coleman-Mcrmin- Wagner theorem, which 
states that in two-dimensional systems the con- 
tinuous chiral symmetry must be manifest for all 



T > 0. Next-to-leading order calculations pL[| 
demonstrated that the the model undergoes a 
Berezinskii-Kosterlitz-Thouless (BKT) transition 
H in accordance with the dimensional reduction 
scenario. In || it was shown that the thermal 
transition of the ^-symmetric GNM3 model be- 
longs to the 2 — d Ising universality class. It is 
easier to visualize the BKT scenario if we use the 
"modulus-phase" parametrization a + in = pe z6 . 
In two spatial dimensions logarithmically diver- 
gent infrared fluctuations do not allow the phase 
9 to take a fixed direction and therefore prevent 
spontaneous symmetry breaking via (8) 7^ 0. The 
critical temperature Tbkt is expected to sepa- 
rate two different chirally symmetric phases: a 
low T phase, which is characterized by power law 
phase correlations (e^^e - * 8 ^) ~ x -71 ^' at dis- 
tances x 3> 1/T but no long range order (i.e. 
a spinwave phase where chiral symmetry is "al- 
most but not quite broken" ) , and a high T phase 
which is characterized by exponentially decaying 
phase correlations with no long range order. In 
other words for < T < Tbkt there is a line 
of critical points characterised by a continuously 
varying < v( T ) < 3- However, since the ampli- 
tude p is neutral under U{1), the dynamics of the 
low temperature phase do not preclude the gen- 
eration of a fermion mass rrif cx p whose value 
may be comparable with the naive prediction of 
the large- Nf approach. In (|] it is shown that in 
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Figure 1. Order parameter |$| vs. (3 for L t — 4 for 
both C/(l) and Z2 symmetric models. 
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Figure 2. Susceptibility % vs - 



0.60 



GNM3 with T > 0, next-to-leading order correc- 
tions cause fermion mass generation to occur for 
T<T„, where T* >Tbkt, i-e. it predicts a "pseu- 
dogap" phase in which the system is non-critical. 
In the next section we present numerical evidence 
in favour of these scenarios || . 

2. Results 

The model was formulated on the lattice using 
the staggered fermion formulation and the simu- 
lations were performed using the standard hybrid 
Monte Carlo algorithm. Details concerning the 
lattice action and the algorithm can be found in 

In order to study the behaviour of the chiral 
symmetry at T > in the absence of a fermion 
bare mass in the Lagrangian the best thing to 
measure is an effective "order parameter" |$| = 
•\/E 2 + LT 2 , which is a projection onto the direc- 
tion of $ Q = (E, II) separately for each configu- 
ration. 

In Fig.|l| we plot |$| versus [3= 1/g 2 for different 
lattice sizes together with the results from simula- 
tions of the model with a Z2 chiral symmetry. At 
the bulk critical coupling (3 b c ulk = 1/g 2 w 0.86 § 
the lattice spacing becomes zero and T— > 00. It 
is clear that the order parameter of the Z2 model 
is independent of the lattice size until just be- 



fore the transition at (3 — f3^ 2 =0.565(3), whereas 
in the [7(1) model |$| has a strong size depen- 
dence for a large range of values of j3, i.e. it de- 
creases rapidly as the spatial volume increases in 
accordance with the expectation that chiral sym- 
metry should be restored for T > 0. The finite 
spatial extent L s provides a cut-off for the di- 
vergent correlation length and according to the 
BKT scenario the slow decay of the correlation 
function (e^e - *^ )) with exponent r;(r)<0.25 
for T < Tbkt ensures a non-zero magnetization 
even in a system with very large size. In Fig.^ 
we plot the susceptibility of the order parameter 
X = V((\&\ 2 ) — (I'&l) 2 ), measured on lattices with 
different spatial size. We observe from this figure 
that: (a) there is a phase transition since the peak 
of x diverges as we increase L s ; (b) the transition 
occurs at a much smaller temperature than the 
critical temperature of the ^-symmetric model, 
because the IR fluctuations are stronger in the 
continuous symmetry case and (c) in the low T 
phase the susceptibility has a a stronger size de- 
pendence than in the high T phase and the errors 
at low T are much larger than at high T. This 
is strong evidence that the system in critical in 
the low T phase in accordance with the BKT sce- 
nario. We also measured the specific heat Cy, 
which we calculated from the fluctuations of the 
bosonic action Sb = \ ^2 x [a 2 {x) + tt 2 (x)]. C v is 
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Figure 3. Fermion screening mass M a vs. /?. The 
horizontal line shows the lattice lowest Matsubara 
mode uo- 



given by C v = fr((Sf) - (St) 2 ). We observed 
that C v has a broad peak at j3 ~ 0.50 and it does 
not show any divergent behaviour or significant 
finite size effects. This is also consistent with the 
BKT scenario, according to which a = —2. 

Finally, we measured the fermion screening 
mass M S (T) from the exponential decay of the 
spatial correlator along one of the spatial lattice 
axes. The results extracted from the simulations 
of the U(l) and models on lattices with differ- 
ent L s , are shown in FigJ|. If we assume that the 
temperature effects are absorbed into a temper- 
ature dependent pole mass M(T) and a "renor- 
malization" of the speed of light coefficient A(T) 
then we have M%(T) = A{T)[M 2 {T) + wg]. We 
can infer from Fig.|3j that the screening mass is 
independent of L s for almost all values of (3. In 
the Zi case we ascribe mass generation to ortho- 
dox chiral symmetry breaking. The similarity of 
Ms (1) and Mf 2 at low T where Ami suggests 
that the pole mass M (T) is nonzero. By compar- 
ing Fig.|| with Fig.|| we infer that the fermions 
are massive in the phase where both chiral sym- 
metry is restored and the system is not critical. 
This implies the existence of a pseudogap phase 
for (3 BKT < < /3f 2 . The deviation of M S {T) 
from u for (3 > 0.56 could be attributted to ei- 
ther A(T) > 1, M (T) > or discretization effects. 



To understand the physics in this regime these ef- 
fects need to be disentangled via a detailed study 
of the fermion dispersion relation closer to the 
continuum limit on lattices with L t >8. For this 
reason we cannot at this stage identify a second 
phase transition where mass generation switches 
off. 

3. Summary and Outlook 

The numerical results presented in this report 
provide strong evidence that the nonzero tem- 
perature phase structure of the t/(l)-symmctric 
GNM3 is consistent with the BKT scenario. We 
have also shown that there is a pseudogap phase 
i.e., the fermion mass remains nonzero above 
Tbkt- We are currently extending our work in 
various directions. We wish (i) to increase our 
statistics in order to extract the exponents T)(T); 

(ii) to study the fermion dispersion relation and 

(iii) to understand the Nf dependence of our var- 
ious results. 
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